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ABSTRACT 

Diffusion-limited reaction A+A— >inert with anisotropic hopping on the d ~ 1 
lattice, is solved exactly for a simultaneous updating, discrete time-step dynamics. 
Diffusion-dominated processes slow down as the anisotropy increases. For large times or 
large anisotropy, one can invoke the appropriate continuum limits. In these limits the 
effects of the anisotropy on variation of particle density can be absorbed in time rescal- 
ing. However, in other regimes, when the discreteness of the time steps is nonnegligible, 
the anisotropy effects are nontrivial, although they are always quite small numerically. 
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1. Introduction 

Recent studies of fluctuation effects in reaction-diffusion systems have been focused 
largely on exact or asymptotic continuum-limit analyses of low-dimensional, d = 1 
and d — 2 models of the simplest reactions such as A+A— i^inert, A+A— i>A, etc. [1- 
12]. Particle input, fragmentation, and multiparticle reactions have been considered as 
well [4,8,12-13]. Fluctuation effects are significant in low dimensions; diffusion-limited 
reactions then show variation of the particle density at large times which is different from 
the predictions of the chemical rate-equation approach (mean- field theory) . In this work 
we initiate studies of the effects of anisotropy in the hopping rates on lattice reaction- 
diffusion models. Specifically, we derive exact results for the reaction A-l-A— >inert, in 
d=l. 

Systems of hard-core particles show interesting novel density fiuctuations associ- 
ated with hopping rate anisotropy, such as shock-waves, etc. Recent studies [14-16] 
have yielded a variety of information pertinent to these phenomena, with most explicit 
results limited to d = 1 models. We note however that the few simplest, solvable 
reaction schemes in d = 1 involve instantaneous reactions such as two-particle annihi- 
lation. These reactions are termed diffusion-limited. The hard core property cannot be 
combined with annihilation on encounter unless the latter occurs with probability < 1. 
The partial-reaction probability models can only be studied numerically [10,13]. Added 
hopping anisotropy may also have a profound effect on reactions in nonuniform systems 
with, e.g., point-like sources of particles, sinks, walls, etc. Investigation of these effects 
is, however, outside the scope of the present study. 

Here we focus on the third interesting effect of the hopping rate anisotropy: slowing 
down of the diffusion. Indeed, let us consider a single particle, hopping on a line with 
the rate H per unit time, and with the fixed step length £ to the right and to the left. 
Let the hopping rate to the right be (1 -|- a)H/2, while that to the left be (1 — a)H/2. 
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Here we assume — 1 < a < 1. In the long-time regime when the continuum description 
can be used, the particle diffusion constant D{a) will be 



Of course the diffusion is superimposed on the drift with velocity a£H. Note however 
that while the time scale of the drift is set by the original time variable t, the time scale 
of the diffusional spread is set by the dimensionless variable 



where the coefficient 4 was introduced for convenience. Thus, diffusion-dominated pro- 
cesses get slowed down as the anisotropy increases. 

Our reaction-diffusion model for the d = 1 reaction A-l-A— s>inert, will be defined 
in Section 2, within the cellular-automaton, simultaneous updating formulation [11], 
with discrete time steps At = 1. However, we will demonstrate (Section 4) that the 
model approaches the continuum-time reaction-diffusion kinetics when the "natural" 
time steps. At = (1 — a^)HAt, vanish in the extreme anisotropy limit \a\ 1. Exact 
solution will be obtained for all a values, in Section 3. The a = case corresponds to the 
isotropic reaction solved in the cellular-automaton variant in [11]. The |a| = 1 limiting 
solution turns out (Section 4) to be identical with the result obtained in the isotropic 
case for the continuum-time Glauber-dynamics variant of the reaction A-l-A— >inert, 
solved in [5,9]. Section 4 also contains a summarizing discussion. 



D{a) = (1 - a'^)e^H/4 = (1 - a^)D{0) . 



(1.1) 



T(a) = AD{a)t/f = (1 - a^)tH , 



(1.2) 



2. Definition of the Model 

We consider particles hopping on a linear x-lattice of spacing £. The dynamics is 
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defined to have the evolution in time steps At = 1. At a given time t = 0, 1,2, . . ., 
we assume that the reacting particles occupy, with certain average density p{t), odd 
lattice sites, x/£ = ±1,±3, at even times {t = 0,2,4,...), and even lattice sites 
{x/£ = 0,±2, . . .) at odd times t = 1,3, . . . . There can be or 1 particles at each 
site (of the relevant sublattice). At each time step t t + 1, each of the particles 
randomly hops one lattice spacing to the right {+£x) or to the left {—£x), with respective 
probabilities (1 + a)/2 and (1 — a)/2. The hopping events are independent. However, 
at those sites which at time t + 1 received 2 particles, the pairs "annihilate," i.e. the 
occupation number is immediately reduced to 0. It is obvious that this rule decouples 
the even and odd sublattices at alternative time steps, — the reason for our restriction 
to the even/odd sublattices as described earlier. 

Quantities such as the particle density per site p{t) are calculated by averaging not 
only over the hopping direction choices but also over the initial particle distribution at 
time t = 0. In order to have an exactly solvable model, here as for the isotropic case 
a = [9,11], one has to consider a particular choice of the initial particle distribution: 
randomly occupied (odd) sites with the average initial density p(0) = 1/2 (per odd- 
sublattice site). 

The root to exact solutions of annihilation models A+A^inert, involves the trans- 
formation [5,9,11] from the hopping particle dynamics to the dual-lattice stochastic spin 
dynamics. In order to introduce the dual lattice systematically, let us describe the orig- 
inal particle system in the following somewhat unusual way: the occupancy variables 
mj(t) will take values —1 or 1, for occupied sites x = j£ or empty sites, respectively. 
Thus the particle density is 



p{t)^{[l-m,{t)]/2) 
-4- 



(2.1) 



Note that the density does not depend on a; = j£ because the initial conditions and 
stochastic hopping rules, — average over both being denoted by (...), — are transla- 
tionally invariant. 

We now introduce the dual-lattice spins <Tj{t) = ±1, on the even sublattice for even 
times, and on the odd sublattice for odd times. Thus, these spins are located in the 
interstices x = j£ between the allowed particle-location sites at a; = {j ± 1)£. Formally, 
the spin values are defined by the infinite products 

oo 

aj{t) = l[mj+2i+i{t) , (2.2) 

i=0 

although with some care in definitions, products up to some fixed mj(t) could be used 
as well. This definition implies 

(t) = {t)<^j+i (t) > (2-3) 

so that the density is given by 

2p(t) = 1 - ((7,_i(t)a,+i(t)) = 1 - G2it) = Goit) - G2{t) , (2.4) 
where we define the two-point spin-spin correlation function 

Gn{t) = {ai{t)ai+n{t)) , (2.5) 

which only depends on the separation \n\ due to translational invariance. We will 
consider this function for n = 0, 2, 4, ... . 

While the stochastic dynamics of the original occupancy variables mi{t) is nonlin- 
ear, it reduces to a linear dynamical rule for the dual spins. This can be shown directly 
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by first introducing the particle dynamics. However, the steps involved are cumber- 
some. A simpler approach is to define the linear spin dynamics first and check that the 
particle hopping and annihilation are correctly represented by it. Thus, let us introduce 
random variables ^i{t) taking on values 1 or 0, with respective probabilities (1 — a)/2 
and (1 + a)/2, where a is the hopping bias. Then the dynamical rule is 

ai{t + 1) = [1 - Ut)] <^i-i{t) + Ut)<^i+i{t) ■ (2-6) 

Having i as a spin site at time t + 1 suggests that at time t the site i was a possible 
hopping-particle location site; recall our sublattice rules. We now identify the random 
variable values ^i{t) = or 1 as the decisions to move the particle to the right or to the 
left, respectively, provided of course that there is a particle at the site x — i£ bX time t. 
Let us also define more natural particle occupancy variables, 

ni{t) = [l-mi{t)]/2, (2.7) 

which take on values (empty site) or 1 (occupied site), and average to p{t). 

When we calculate mi(t + 1) in terms of the ruj variables at the preceding time, 
t, by substituting (2.6) in the t + 1 version of (2.3), the relation is no longer linear. 
The source of the nonlinearity is the term proportional to (7i-2(7i-^2 = fni-inii+i. After 
some algebra, we get, in terms of the variables, 

riiit + 1) = ^i+i{t)ni+i{t) + [1 - ei-i(t)] m-iit) - 2^i+i{t)ni+i{t) [1 - ei-i(t)] ni_i(t) . 

(2.8) 

It is obvious that this nonlinear dynamics corresponds to the particle-hopping rules as 
specified originally. Indeed, the first term adds 1 to nj(t -|- 1) if there was a particle 
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at site i + 1 at t, and if it hopped to the left. The second term adds 1 if there was a 
particle at i — 1 and it hopped to the right. Finally, the last, "reaction" term, which 
is (—2) times the product of the first two, "diffusion" terms, sets the resulting value to 
zero if both diffusion terms were 1, i.e., if two particles hopped to i in the time step 
t — > t + 1. Unfortunately, we found no other nontrivial quadratic dynamical rules that 
can be linearized by the dual-spin transformation or its simple generalizations. Thus, 
the reaction A+A— >inert is the only one solvable by the dual-spin method. For instance, 
if we replace (—2) by (—1) in (2.8) to have the reaction scheme A-l-A— >A, then the dual- 
spin method fails. The latter reaction can be solved by another method, however; see 
[12]. 

Finally, we note that the hopping rate of a single particle is = 1. Thus, the 
dimensionless time variable r is given by 

T={l-a^)t. (2.9) 

3. Exact Result for the Particle Density 

The initial particle distribution was selected to correspond to the random ±1 as- 
signment of the initial {t = 0) spin values. Averaging (2.6) over the stochastic variables 
^ then leads to the conclusion that (c"j(t)) = for all times t > 0. The two-point corre- 
lation function, however, is nontrivial. Forming the product suggested by the definition 
(2.5), we obtain, after averaging and some algebra, 

Gnit + 1) - Gr^it) = [Gn-2{t) - 2Gn{t) + Gn+2(t)] • (3.1) 

This relation applies only for n = 2, 4, 6, . . ., and it must be solved subject to the 
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boundary condition 



Go{t > 0) = 1 , (3.2) 

as well as the initial condition 



G'n>o(0) = 0. (3.3) 
Let us introduce the generating functions 

oo 

Bn{v) = J2^'Gn{t), (3.4) 

as well as the generating function for the density, 

oo ^ 

Riv) = Yl "V(i) = 2 [^o(^) - B^iv)] , (3.5) 
where the last equality follows from (2.4). Note that we have 

Bo{v) = , (3.6) 

1 — V 

while the n > generating functions satisfy the recursion 

B 1 — 0? 

Bn= {Bn-2 - 2Bn + Bn+2) ■ (3.7) 

V 4 

Each term in (3.7) follows directly from summation of the appropriate term in (3.1). 

This second-order difference equation is quite standard. Only one out of two inde- 
pendent solutions behaves regularly for large n, in the vicinity of v = 0. The solution 
takes the form 
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Bn{v)^Bo{v) [b{v) 



i/2 



(3.8) 



where the appropriate function b{v) is found directly from (3.7), 



(1 — a^) V 



1- W1 + 



[l-a^)v 
1-v 



The generating function for the density then foUows as 



(3.9) 



R{v) = 



1 



\/l — a^v 



1 



(3.10) 



where we made several algebraic steps to prepare the resulting expression for the Taylor 
series expansion to yield the density: 



(1 - a') Pit) 



{2t + 2) 



(2n)!(2t-2n)!a2"+2 



l)!f "^o2^'+'^!(n+l 



(n+l)![(t-n)!] 



(3.11) 



4. Special Limits and Discussion 

In this section we consider three limiting cases in which the solution (3.11) assumes 
less cumbersome forms. First, we note that for the symmetric- hopping case, a = 0, 
only the first term in (3.11) survives. Since r = t for a = 0, we have 

P,=o(t) = + , . (4.1) 

^ ' 22^+2 [(r + l)!]' ^ ' 

This form provides a convenient analytic continuation to all r > 0, shown as the lower 
solid curve in Figure 1. However, the reader should note that the actual time-dependence 
for |a| 1 is always in discrete time steps. 
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Next, consider the limit \a\ — > 1, i.e., extreme anisotropy and slow diffusion. The 
"natural" time steps At = (1 — o^) then become small and the time-dependence ap- 
proaches a continuous function of r. As a result, the discrete-time difference on the 
left-hand side of (3.1) can be replaced by the derivative (1 — a^)dGn{T)/dT. The re- 
sulting equations are then identical, up to minor notational changes, to those obtained 
in the continuum-time Glauber-dynamics formulation [5,9] of the (isotropic) reaction 
A-l-A— >inert. They can be solved, for instance, by the Laplace Transform method. We 
only present the final expression for the density. 



P|a|=i(r) = le-^/' [/o(r/2) + /i(r/2)] , (4.2) 



where /„ are the standard Bessel functions. The upper solid curve in Figure 1 corre- 
sponds to (4.2). 

Finally, we note that for large times both the temporal and spatial variation of 
the two-point correlation function can be regarded as slow and differences replaced by 
derivatives. Specifically, the second-difference expression on the right-hand side of (3.1) 
is replaced by the second spatial-coordinate derivative. The detailed derivation of this 
continuum limit and solution for two-point correlations are presented, e.g., in [11]. In 
this limit, the hopping anisotropy is fully absorbed in the rescaling of time from t to t. 
We only quote the result for the density 



p(t » 1) ~ . (4.3) 

VTTT 



This function is shown by the dashed line in Figure 1. It is asymptotic to the solid 
curves for large r, as well as to the full result (3.11). 

The full expression (3.11) is more complicated than the limiting cases just discussed. 



The density, defined for discrete r-values, — integral multiples of (1— a^), — lies between 
the two limiting solid curves shown in Figures 1 and 2. For fixed t, the r values decrease 
as \a\ 1. Thus, the crossover from the a = discrete values to the |a| = 1 continuous 
T-dependence occurs for fixed t as shown in Figure 2. The points in Figure 2 correspond 
to t = 2 and 

a = 0.00, 0.05, 0.10,..., 0.95 ; (4.4) 

note that the symbols for the points corresponding to a = and 1/20 largely overlap in 
the figure. 

In summary, we derived exact results for a two-particle annihilation reaction with 
anisotropic hopping on a d = 1 line. The main effect of the anisotropy is to slow down 
the diffusion-dominated processes, and it can be largely absorbed in the redefinition 
of the time scale. However, in the regimes where the discreteness of time steps is 
nonnegligible, i.e., for short times and nearly isotropic hopping, the functional form of 
the density depends on the hopping anisotropy parameter a explicitly. The discrete- 
time-step reaction proceeds slightly faster than the continuous-time reaction, though 
the difference is numerically quite small. 

The author wishes to thank M. Barma, M.L. Glasser and S. Nicolis for helpful 
comments and suggestions. 
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FIGURE CAPTIONS 



Figure 1: Particle density as a function of the dimensionless time variable r. Lower 
solid curve: a = 0, the actual densities correspond to integer r values; see Section 4. 
Upper solid curve: a = ±1. Dashed curve: the r » 1 asymptotic form. 

Figure 2: The points illustrate the crossover from a = to |a| = 1, for a fixed time, 
t = 2. The a values are given in Section 4; see (4.4). The solid curves are the same as 
in Figure 1. 
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